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Abstract 

A bi-layer quantum frustrated antiferromagnet is studied using an effective 
action approach. The action derived from the microscopical Hamiltonian has 
the form of the 0(3) non- linear sigma model. It is solved in the mean field 
approximation with the ultraviolet cut off chosen to fit numerical results. 
The obtained phase diagram displays a decrease in the critical value of inter- 
layer coupling with increase of in-plane frustration. The critical point for a 
single-layer frustrated antiferromagnet (the J\-Ji model) is estimated to be 
J 2c = 0.19 Ji . 
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In the past few years, a lot of attention has been given to order-disorder transitions 
in two-dimensional antiferromagnets. Much of the interest stems from relevance of the 
problem to high-temperature superconductors where such transitions occur. Although in 
superconductors the transition is driven by mobile holes, it can also be studied in a frustrated 
antiferromagnet. The antiferromagnet with frustrating interaction between next nearest 
neighbours (the J1-J2 model) is the simplest case and has often been studied.! 

In this paper, I study a bi-layer frustrated antiferromagnet. The bi-layer problem has 
emerged from the studies of YBa 2 Cu 3 6+ ,j; which consists of pairs of close Cu0 2 planes. Ex- 
perimentally, there are a number of differences between bi-layer and single layer compounds.i 
Recent infrared measurement si suggest that the antiferromagnet interlayer coupling within 
bi-layers may be as strong as 60 meV (about half of the in-plane coupling) in which case it 
must be taken into account in any theoretical study. Several authors have discussed pair- 
ing in bi-layer models.^ This work is restricted to studying magnetic order in the static 
ant if erromagnet . 

The Hamiltonian is taken in the form 



H — J\ XX^lr ' Sir' + S 2r • S 2r ') + J2 XX^ lr ' Sir' + S 2r • S 2r ') (1) 

NN 2N 

Ir ■ Sir' + S 2r ■ S 2r ' Sir ' S 2r , 



3N 



where the summation is performed over a square lattice, NN denotes nearest neighbours, 
2N (3N) stands for second (third) neighbours, and the two planes within the bi-layer are 
labelled as 1 and 2 (it is assumed that all interactions J are positive). 

In case J± = the Hamiltonian ([!]) reduces to a frustrated Heisenberg antiferromagnet. 
Increasing J 2 and J3 leads to a quantum phase transition into a disordered phase. The 
numerical value of the critical coupling J 2c (at J 3 = 0) is, however, known very poorly. 

The case of J 2 = J3 = (a bi-layer antiferromagnet) has also received a lot of 
attention.Ull Availability of reliable numerical resultslE for this case makes it a conve- 
nient mo del§ for testing theoretical results. The bi-layer antiferromagnet has also been 
used! to explain some experimental results in bi-layer superconductors. However, the value 
of interlayer coupling required for this (J± > 2.5Ji) is too large to be realistic. Apparently, 
both in-plane and inter-plane effects are important. In the presence of in-plane frustration 
the critical value of interlayer coupling may be drastically reduced. 

In both the bi-layer antiferromagnet and single layer frustrated antiferromagnet, the 
mean field spin wave and mean field Schwinger boson theories significantly overestimate the 
region of stability of the Neel ordered phase. The main reason for this appears to be that 
these theories assume long range order from the start and do not include some types of 
fluctuations (see discussion of longitudinal spin fluctuations in bi-layer antiferromagnet in 
Ref. 0)- 

In this study, I use an effective action approach which is based on mapping the micro- 
scopic Hamiltonian onto the quantum 0(3) non-linear sigma model and which has been 
extensively applied to two-dimensional antiferromagnets. Using a 1/N expansion of the 
non-linear sigma model, Chubukov, Sachdev, and Jinwu Ye0 have recently presented a 
comprehensive analysis of the general properties of clean two-dimensional antiferromagnets 
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in the vicinity of the order-disorder transition. Good agreement with numerical and exper- 
imental data has been obtained. 

In this work, I first derive the effective action for the bi-layer frustrated antiferromagnet. 
It is observed then that for the purposes of the study, the action can be reduced to that of 
the standard one-band non-linear sigma model. The phase diagram is then found using the 
mean field (saddle point) solution with the ultraviolet cut off (or, equivalently, the critical 
coupling) adjusted so as to fit numerical results for the bi-layer antiferromagnet. 

There are several ways to obtain the action. I follow the derivation based on coherent 
state representation.0 Each spin is represented by a single unit vector N. Then all spin 
configurations are expressed as a combination of four fields (two fields per layer) 



N lr 
N 2r = 



77rni r (l 

-r] r n 2r (l 



2 LL) 1/2 + «L lr , 



(2) 



where 



n7 



n 



1, ni ■ Li = n 2 ■ L 2 = 0, a 2 L 2 2 < 1 



r] r = ±1 for the two sublattices, and a is the lattice constant. The unit length fields rij r in 
the classical limit describe orientation of the local Neel ordering and the fields Lj r represent 
small local fluctuations. 

Now the Hamiltonian is expressed using Eq. @ and the continuum limit is taken by 
making an expansion in gradients of ni and n 2 and powers of Li and L 2 



H = J rf 2 r[ip sl [(V ni ) 2 + (Vn 2 ) 2 ] + S'J^n, - n 2 f 
+4S 2 (J 1 + U±){Ll + L 2 ) + S 2 J ± U • L 2 



(3) 



where p s \ = S 2 (Ji — 2J 2 — 4J 3 ) is the bare spin stiffness. 

Decomposition of the partition function is a straightforward extension of the procedure 
for the single layer antiferromagnet0 and leads to 

Z = J Pn 1 I?n 2 I?L 1 I?L 2 5(n 2 - l)5(n 2 - 1) exp(S n ) (4) 

with the action 

S n = S B - J d 2 rj^dr[H(T) + iSU ■ (m x cUx) + iSL 2 ■ (n 2 x <9 T n 2 )], 

where Sb is the residual Berry phase which is zero for all smooth spin configurations and 
which is ignored hereafter. Finally, the fields Li and L 2 are integrated out of Eq. (U) 

Z = J Vn x Vn 2 8{n{ - l)5(n 2 - 1) 

ip sl [(V ni ) 2 + (Vn 2 ) 2 ] + 5 2 J ± ( ni - n 2 ) 2 



x exp< 



c/3 



dr 



+ l(Ji + jJ±) l { n i x + n 2 x <9n 2 ) 2 + ij x ^ni x 9ni - n 2 x dn 2 ) 



(5) 
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At J± = 0, the action (|5]) represents two independent sigma models each possessing a 
Goldstone mode at T = due to spontaneous symmetry breaking. At finite J±, one mode 
acquires a gap A oc yJ±J\. There is no general solution in this case because of the presence 
of the term (ni x dux) ■ (ri2 x <9ri2) 2 ) in the action. However, at larger J± and T <C y/J±Ji 
the action can be simplified due to the fact that it is dominated by configurations with 
ni « 112. In this regime of coupled planes, I restrict myself to the Goldstone mode and set 
n = ni = ri2 in Eq. (^). The action takes the form 

Z = J Vn5{n 2 - 1) exp j ~ J d 2 r jf" dr [( Vn) 2 + (d T n) 2 ] j , (6) 

where time has been rescaled (cr — ► r ) and the coupling constant is g — c/p s with the spin 
wave velocity 

c = 2S[2(J 1 - 2J 2 - 4J 3 )(Ji + \ J ± )} 1/2 (7) 
and the bare spin stiffness 

p s = 2S 2 (J 1 -2J 2 -4J 3 ). (8) 

From Eq. (0) for the bi-layer antiferromagnet the spin wave velocity at the critical point 
J± = 2.5Ji is c c = I.8J1. It is in good agreement with numerical resultsS which suggests 
that 1/S corrections are negligible at this point. The spin stiffness Eq. (§) is twice as large as 
in the regime of independent planes. Thus when J± is introduced the system first becomes 
more "classical" (in other words, the effective spin S* increases) and then evolves towards 
the quantum disordered phase. This effect was observed in earlier studiesiffi at a fairly 
small value of J± ~ O.lJi. With in-plane frustration present the crossover will occur at even 
smaller values of J±. (In a recent exact diagonalisation study of the bi-layer t-J model,El it 
was noticed that the two layers become essentially correlated at w 0.2 J.) 

Equation (||) represents the quantum 0(3) non-linear sigma model. Not being interested 
in details of the critical behaviour, I use the simplest mean field solution which is exact for 
the O(oo) non-linear sigma model. An ultraviolet cut off must be introduced in momentum 
integration and a Pauli-Villars cut off is a convenient choice.^ In the mean field approxima- 
tion, there are spin-1 excitations with a gapped spectrum and no damping. Straightforward 
calculations show that the gap m is determined form 

(9) 

where g c = 4ir/A is the critical coupling (A is the cut off). The system is ordered at T = 
when g < g c . Knowing the critical point at J 2 = J3 = allows to eliminate uncertainty from 
the cut off. Using Jj_ c = 2.5J, the cut off is found to be A = l.ln. It is easy to see that the 
critical line g = g c has the form 

J hJ 1 

j_ 2 J 2 -4J 3 = j,? T - (10) 

J + l J .Lc 











( sinh — ] 
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V 2TJ 
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TABLE I. The critical coupling in the J\-Ji model obtained by different methods. 



method 

Green functions 3, 
renormalization group b 
mean field bond operators 
series expansion d 
present work 
exact diagonalisations 

a Ref. |e 
b Ref. [17 
c Ref. p 
d Ref. [IJ 
e Ref. Eo 



J 2c 

0.12 
0.15 
0.19 
0.33 ±0.08 
0.19 
> 0.34 



While Eq. (Q) is only the mean field solution and is not quite accurate, Eq. ([U]) is independent 
of 1/N corrections. The phase diagram is presented in Figs. 2 and 3. Figure 2 is the phase 
diagram for J 3 = 0. In the absence of interlayer coupling, the critical value of J 2 is found 
to be J 2c = 0.19Ji. It is compatible with the results obtained by other methods except 
for numerical results (see Table I). However, numerical studies^ face severe size restrictions 
and are not very reliable quantitatively. Figure 3 is the phase diagram for J 2 = 0. Both 
Figs. 2 and 3 show that in the presence of an interlayer coupling smaller values of in-plane 
frustration are sufficient for the disordering transition. 

In conclusion, I have studied the bi-layer frustrated antiferromagnet using the effective 
action approach. It was demonstrated that the essential physics can be described by a one- 
band non-linear sigma model. The region of stability of the Neel ordered phase has been 
identified The critical value of next nearest neighbour interaction for the J1-J2 model is 
estimated to be J 2 = 0.19Ji. 
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FIG. 1. A scheme of the bi-layer antiferromagnet. 



FIG. 2. The phase diagram for J3 = 0. 



FIG. 3. The phase diagram for J2 = 0. 



6 



